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Abstract. An Af'-particle system with stochastic interactions is con- 

**< ■ sidered. Interactions are driven by a Brownian noise term and total en- 

_ ■ 

ergy conservation is imposed. The evolution of the system, in veloc- 
ity space, is a diffusion on a (3N' — l)-dimensional sphere with radius 
i-G ' fixed by the total energy. In the N' — >• °° limit, a finite number of veloc- 

ity components are shown to evolve independently and according to an 

i r*j Ornstein-Uhlenbeck process. 
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1.. Introduction 

The ensemble average of the stochastic evolution of a Brownian test par- 
ticle immersed in a heat bath of constant temperature T is given by the 
Fokker-Planck equation with prescribed coefficients of diffusion and linear 
friction HUO30I . The solution to this equation is well known and gives the 
evolution of the ensemble's probability density function in velocity space. 

Kies sling and Lancellotti [KL06] have shown that, if the coefficients 
of this partial differential equation for the Ornstein-Uhlenbeck process are 
chosen to be constant functionals of the solution itself, the equation can be 
reinterpreted to describe the kinetic evolution of an isolated A^-particle sys- 
tem with stochastic interactions satisfying total mass, energy and momen- 
tum conservation. In their work, these authors consider an infinite ensemble 
of independent and identically distributed (i.i.d.) random iV'-dimensional 
vectors {V a }^ =1 , each representing a possible microstate of the system in 
velocity space (in M? xN ), while particle positions are considered to be uni- 
formly distributed over a periodic box. By a suitable transformation on the 
V a 's, they show that the analysis of this A^-particle system with energy and 
momentum conservation can be done in terms of an A^'-particle system with 
only energy conservation and N" = N' — 1. 

Our purpose in this work is to analyse the evolution of the same parti- 
cle system with stochastic interactions driven by this Brownian noise term, 
but from a many-body perspective, i.e. we look at a single realization of 
the system and study its evolution pathwise. According to the arguments 
of HKL06H . it suffices to study the system with only the energy conserva- 
tion requirement, and so we do here, considering the "velocities" of this 
effective model as if they related to actual particles (in the spirit of Kac' 
one-dimensional historical model [Kac56J). This model also belongs in the 
class of Kac systems in the terminology of HCCL03H . 

In a 3-dimensional setting, energy conservation bounds the evolution of 
the system, in velocity space, to an (N — 1) -dimensional sphere with ra- 
dius defined by initial energy (with N = 3N"). Therefore, the evolution of 
the system in velocity space can be modeled by the Brownian motion on 
the given sphere. Thus, we take advantage of the works of Stroock [StTTJ 
and Brillinger HBr971 , where the random evolution of a single particle on a 
sphere was modeled with stochastic integrals, to write the proper differential 
equations for the evolution of the A^-dimensional vector V of all the velocity 
components for the N' particles. Singling out one component of V (call it 
Vi), we show it to converge, in the N — > °° limit, to an Ornstein-Uhlenbeck 
process HUO30L which is commonly used to describe the evolution of a 
single particle under "white noise" and friction. 
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2.. Outline and strategy 

We look at the diffusion of a 3 -dimensional A^'-particle system which 
obeys conservation of energy. This requirement bounds the evolution of the 
system, in velocity space, to a (3N" — 1) -dimensional sphere with radius 
defined by initial energy. The dynamics of the system is given by a single 
stochastic differential equation for a (N = SA^O-dimensional velocity vector 
(V) driven by an A^-dimensional Brownian noise term introduced in Sec. [3T] 

A quick look at this evolution equation shows that a single component 
of V evolves independently of the remaining directions according to an 
Ornstein-Uhlenbeck process driven by a single noise term along the same 
direction, when this component is small enough. Our focus, however, is on 
studying the limiting process for the components of V when these are of 
order of unity. Formally, we define a reference one-dimensional Ornstein- 
Uhlenbeck process, and show that the evolution of a single component of 
V converges (in probability) to this reference process as N — > °°. This is 
properly stated in Sec. |4]It is proved in Sec. |5] 

This proof is easily extended to any finite ^-dimensional process [n < N) 
constructed from n components of the original diffusion. Thus, we show 
that these processes are independent and identically distributed (i.i.d.) in 
the N — > °° limit, thereby proving the propagation of initial independence of 
n particles, viz. the "propagation of molecular chaos" or of "Boltzmann's 
property" ( [|Kac56[ |Kac59IL see the recent review by Mischler and Mouhot 
[MM1 1J for an extensive discussion), in Sec. [6] 

In the appendix, we sketch the answer to a question raised by Kies sling 
and Lancellotti [KL06] on the explicit characterization of the many -body 
stochastic process. 

3.. Diffusion on a sphere 

Consider a system of N" particles. The particles exchange momentum 
according to some interaction between them, and total energy is assumed to 
be conserved. 

We define the ^"-dimensional vector of interest 

(i) \" = (vi',vl...,v^), 

where V" represents the velocity of the z-th particle in M?. Conserva- 
tion of energy requires the vector V /; to assume values in the (3N" — 1)- 
dimensional manifold of total energy N"eo 

(2) M™"- l = {v> : £^=AT% 




viz. the (3N" — 1) -dimensional sphere with radius y/2N"eQ. 
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Let N = 3N" denote the number of degrees of freedom of the system. 
Define the ^-dimensional vector 

(3) Y=(V h V 2 ,...,V N ), 

where the V,'s are all the N = 3N" components of the three-dimensional 
vectors Vf. Starting from an arbitrary initial data, the evolution of vector V 
is modeled by the Brownian motion on the sphere, given by the stochastic 
differential equation HBr971 

(4) dV = Aa(V)-odB, 
where A is the noise amplitude, 

(5) ^(V) = %-p 

is the projector onto the direction orthogonal to V, 

(6) B = (B U B 2 ,...,B N ) 

is the standard A^-dimensional Brownian motion, and o d denotes the use of 
Stratonovich calculus [0klO|. Although the Stratonovich view-point seems 
physically natural in the sense that it leads to the same chain rules as classi- 
cal calculus [WZ65], we shall change view -points and prefer working with 
the Ito representation. This representation has some advantages to our pur- 



poses, for instance, Ito integrals are well known to be martingales [|0klO| , 
which is crucial for stochastic analysis and computing expectations. The 
counterpart to this advantage is that Ito integrals do not transform so nicely 
under transformations, e.g. the chain rule for Ito differentials involves sec- 
ond order terms f|0kiol IKPS03H . 



Therefore, we work with the spherical Brownian motion as the solution 
of the Ito differential equation (lSt7Tl[Pr05l 

(7) d\ = Xo(Y)-dB-X 2 ^-^-b(\)dt, 
where 

(8) b(\) 



IV 



2* 



The process © remains on the surface of the sphere, and the last term on 
© may be thought of as pulling the process back onto the sphere HBr97t 
Introducing := 2eo/3, we rescale time to t' = X 2 t/6, the Wiener process 
to B ; = A0~ ! ' 2 B and velocity to V' = _1 ' 2 V, so that, dropping primes, 
the system CO)-® reduces to 

(9) dV = a(V)-dB-— ^-Vdf. 
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4.. Convergence to the Ornstein-Uhlenbeck process 

We single out one component of V to study its evolution. Evidently, this 
evolution will have a contribution from all the other directions (the N — I 
remaining dimensions), so we define 

(10) U = V-Viei, 

where e\ is the unit vector along the direction of the ?'-th dimension. A 
simple rewrite of ® leads to the system 

N-\ 

y 



(11) dVi = <7i(Vi)dBi-— — *i(Vi)d*-D(U,Vi).dBCT, 



(12) dU = ofr(U)-dBtf-^— MU)dt-D(U,Vi)dfli, 



given 



*i(Vi)= 1-iSr, *i(Vi) = i^r, D(D,Vi) = JiU, 



(13) (au(V))ij= 5ij-^ b f/ (U) = r ^, Bu=(B 2 ,...,B N ) 

with the constraint |V| = v^V- 

One obtains a simple estimate of the evolution of Vi by first looking at 
the limit of small values for V\ (Vi ~ 0). In such a limit, 

(14) <Ti(Vi)~l, D(U,Vi)~0. 
So, the resulting evolution equation would be 

(15) dVi-dBi-^^lvidt, 
which reads, in the N — > °° limit, 

(16) dVi-dfii-yd/ 1 , 

which has for solution the standard Ornstein-Uhlenbeck process. 

Let us now look at the general evolution of V\ in the N — V °° limit, i.e. 
we want to know what is the limiting process defined by (fTTT ) when V\ is of 
order one. From the result (|T6l , we can expect this process to, at least, be 
similar to the Ornstein-Uhlenbeck process. Our main result is expressed in 
the following 

Theorem 4..1. LetV\ be a 1 -dimensional component of "V defined by equa- 
tion ([9]). For |V| = y/N, the process V\ with initial data c\ obeying 4771) 
converges, in the N — > °° limit, to the I -dimensional Ornstein-Uhlenbeck 
process with the same initial data and control B\, and convergence is in 
probability in C( [0, T] , M) for any time T > 0. 

Of course, given c\ £ R, the finite-iV model makes sense only for N >c\. 
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5.. Proof of Theorem @TT| 

We start our analysis by defining the reference (or target for the limit) 
1 -dimensional Ornstein-Uhlenbeck process 

(17) dV[ = adB l -^V[dt, a,j8eR. 

coupled to (fTTI) by its involving the same control B\ . With this, we introduce 

(18) v = V l -V{. 

We compute the expectation of the square of this new processes to grasp 
the evolution of Vi compared to the evolution of the reference process. So, 
we start by writing explicitly 
(19) 



v(0=vo+ 



[a 1 (V l )-a)dB l 



N-l 



bi(Vi)-pv{) df-D(U,Vi)-dB 



u 



where vo = v(t = 0). 

We set the parameters a = 1 and /3 = 1/2. Thus, the expectation reads 
(20) 



E[|v(0l 2 l=E 



vo 



* fv(s) Vi(s) 
2 2N 



ds 



Vt 



'Vi. 



o N l Jo N 



which is readily bounded as 
(21) 

3^ffl<E[|vol 2 ] + E 



f fv(s) Vi(s) 



IN 



ds 



+ E 



o N 



V-dB 



Using the property d(B i: Bj) = 5ydf, where (•,•} denotes the cross- 
variation process [0k lO| , and the fact that E [fQf(s)dBj(s)~\ = 0, for any 
function of time /, the last term in the r.h.s. can be rewritten as 

2" 



N 



(22) £E 



(=l 



o N 



VidBt 



N 



(=1 



, -^Vfds 

o N 2 ' 



E 






ds 



where we used Ito isometry in the intermediate step and recalled our con- 
straint, | V| = y/N, in the last step. 
Going back to (|2TI) . we further bound 

E[|v(0l 2 ] 



< E[|v | 2 ] + -E 



V^ds 



2 Jo 



< E[\v \ 2 ] + - E[v(s) 2 ]ds + 



+ l E 



2N 2 4 Nj Jo 



' Vi(s) 

o N 

1\ rt 



ds 



+ E 



,V -IJ 

o N 



E[vl(s)]ds, 
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where we used Schwartz inequality on the second and third terms in the 
r.h.s. To get a closed estimate, we must now control how fast E [V 2 ^)] 
grows in time. For this purpose, a standard Ito calculation leads to 

/ T/2\ 



(23) 


dVf = 2VidVi+[l--^ldf, 


which gives 




(24) 


dE[ ^(,- E ^) d ,<d,, 


so that 




(25) 


t-t rt f2 

/ E\Vi(s)]ds< / (c 2 l +s)ds = cjt + - 



2 
Thus, we find a bound on E[|v(?) | 2 ] independent of the component V\, 

E[|v(0| 2 ] < 3E[|vo| 2 ] + |/ o f E[v(,) 2 ]d, 



'0 

(26) +3 i^4)(* + y: 



1 T IT /"' 

(27) < 3E[|v | 2 ] + — (1 + — )(2c 2 r + r 2 ) + T y o E[v(,) 2 ]d, 

over the interval [0, T] for any T > 0. Gronwall's lemma [Di80| then yields 
the bound 



(28) E[|v(OI1< 



3E[|v | 2 ] + ^(l + ^)(2c 2 r + r 2 ) 



JTt/2 



For any finite t, we may set T = t in this estimate. Therefore, in the 
N — 7- °° limit, we have for finite times 

(29) E[|v(0| 2 ]<3E[|vo| 2 ]e 3/2 / 2 , 

ensuring that, given initial conditions such that vo = 0, the process V\ con- 
verges in probability to the reference Ornstein-Uhlenbeck process in one 
dimension over [0, T] V T < °°, thus proving the theorem. □ 

6.. Motion of n particles 

This result can be generalized if we consider, instead of the single compo- 
nent V\, a finite rc-dimensional vector Ui made up of the first n components 
of V. Define 

n N 

(30) Ui = £y^, u 2 = £ v t iu 

i=\ i=n+\ 
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for a fixed n, and rewrite © as 

(31) dU 1 =a 1 (U 1 )-dB 1 -^b 1 (U 1 )dr-D(U 1 ,U 2 )-dB 2 , 

(32) dU 2 = a 2 (U 2 ) • dB 2 -^lb 2 (U 2 )dr-D r (Ui,U 2 ) ■ dB 1; 
where T denotes the transpose and 



(i%(u,)«. um ,_ Ui 

IVI 



(33) (Oi) mn = 8 mn '" 2 ' - ; b/(U/) = — |j; /=1,2 



(34) (D), 



(Ui) w (U 2 ) n 



|V| 2 

(35) Bi = (5i, ... ,5„), B 2 = (5 n+ i, ... ,B N ). 

From this point of view, we derive the following 

Theorem 6..1. Let \J\ be the n- dimensional component of V evolving 
according to equation 071) . For |V| = y/N, the process Ui vvzYft initial 
data Ui (0) = Ci et° converges, in the N — > °° //m/?, to ?/?e n-dimensional 
Ornstein-Uhlenbeck process with initial data Ci anJ control Bi, anJ con- 
vergence w w probability in C( [0, 7] , IR") /or anv rime T > 0. 

Sketch of proof : The proof parallels the previous one, we need only define 
the reference n-dimensional Ornstein-Uhlenbeck process Uj. For this case, 
the form of the final bound on E[|u(f)| 2 ], with u(t) = Ui(t) - U[(t), is 
identical to (|28l) . replacing v with u and c\ with Ci. 

Note that |ci| is independent of A^ (typically 0(n)), which still leads to 
a vanishing contribution from the second term in the r.h.s. of the equation 
corresponding to d28T >. □ 

7.. Conclusions 

A 3-dimensional A^'-particle system, in which particles interact stochas- 
tically due to a driving Brownian noise term, satisfying total energy (N"eo) 
conservation is modeled by a diffusion process on the sphere M^ . In 
the JV" —7- oo limit, a single velocity component evolves independently of all 
the remaining directions for finite times. Furthermore, this component con- 
verges to a 1-dimensional Ornstein-Uhlenbeck process driven by the single 
noise component coupled directly to it. The role of friction in the usual 
Ornstein-Uhlenbeck model to keep the velocity component bounded is here 
played by total energy conservation and by the curvature of the sphere 
(equatorial bands have a larger area than polar caps). 

Technically, convergence to the Ornstein-Uhlenbeck process is proven 
using standard arguments (see section 5.2 of [0klO|). In the language of 
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IIKPS03L we obtain mean-square convergence for the velocity process in 

An immediate extension is given for the n-dimensional velocity pro- 
cess. In the N — > oo limit, these are shown to be i.i.d., converging to the 
^-dimensional Ornstein-Uhlenbeck process. If at initial time the n tagged 
components have independent data, the subsequent evolution of these n 
components preserves their independence, which amounts to propagating 
molecular chaos. 
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Appendix : Two representations of the diffusion on the 

SPHERE, AND A BINARY INTERACTION INTERPRETATION 

Our representation of the diffusion on the sphere, with the Ito differential 
equation ©, describes the evolution of the velocity vector in terms of N 
independent brownian controls Bj, 1 < j <N, each directly associated with 
one component of V. Kiessling and Lancellotti [KL06] recall the represen- 
tation of the Laplace-Beltrami operator on the unit sphere E> N ~ 1 C M. N in the 
form 

(36) Agw-i = £ (v k d vt -vid Vk f 

\<k<l<N 

and interpret the Fokker-Planck equation in terms of a model for their orig- 
inal N' = \ +N/3 particles with two-body and one-body interactions pre- 
serving momentum and energy. Here we translate this interpretation to a 
simple stochastic process, limiting ourselves to the N/3 "effective parti- 
cles" with interactions preserving only energy. 

Diffusion on the sphere may be viewed as a succession of independent 
infinitesimal rotations. So, consider a family of N(N — l)/2 independent 
processes &ki(-), I < k < I < N, taken to be martingales with initial value 
and cross-variation (Q i y,Q y t/)(0 = 5#5« t/N. Complement Q. to an an- 
tisymmetric matrix, with £2y = — Q,ji, so that the differentials d£lij act as 
infinitesimal rotation generators, and consider the process V' defined by an 
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initial data V'(0) with | V'(0) | = a//V and the Stratonovich differential equa- 
tion 



(37) dV; = Y,VjOd£li 



Its Fokker-Planck operator, acting on measures / on the sphere S^ 1 

with radius a//V, is the sum L = — 5 Li<;<;<A^/Vfti> where Rtj is the vector 
field associated with Qy, and R*: is its adjoint operator on R . By (1371) . the 
vector field is R tj = AT 1 / 2 ^^ - vjd/.). Then R* tj = -N'^id^v'j •) - 

d v ' ( v 'i ') ) = — -fy/> so mat ^ = M^s^-' • A s me generator determines the law 
of a diffusion, this proves that the process defined by (l37l is the same as the 
one defined in Sec. [3] 

Now, comparing directly the differential equations is also interesting. So, 
from the Stratonovich form (l37T i, we deduce the Ito equation for V', 



(38) dv; = £y;dn j7 +i£d(y;,%) 

j i 

= EYJdOy + iEEvMfybfly) 

j z J &J 

(39) = ^yjd^-k^i-^ 1 )^. 



For the first equality, we used the relation between Ito and Stratonovich 
integrals (see e.g. Sec. V.5 in [Pr05]) ; then we substitute the first term of 
the r.h.s. of (1381) into the cross-variation term, and finally we use the explicit 
cross-variation of £1. 

In the final expression (l39l , the drift is exactly the one in (|9]). The first 
term appears as coupling the z'-th component of vector V ; with all other com- 
ponents, and it is linear with respect to V ; while © involves a projection 
matrix o orthogonal to V. As this first term is the differential of a martin- 
gale, say dM'j = ^Vydf^, we now show that it is equivalent (in law) to 
the first term in ® . Indeed, this martingale M' is further characterized by 
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its cross-variation process, which has the differential 

k I 

k I 

(40) = ir^Sujyf-viVj)* 

k 

where one readily recognizes the projector entry 0/,-(V'). 

For comparison, the first term in © defines a martingale M with dM ; = 
£y <Jjj(Y) dBj, so that its cross-variation satisfies 

(41) d(Mi,Mj) = E^OflfcOy/d^fi/) = <7y(V) df 

as follows from the cross-variation of B and the fact that o 2 = o. 

As a result, M' and M have the same cross-variation when generated from 
the same trajectory (assuming M(0) = M'(0) = for definiteness), as befits 
processes V' and V having the same law (incidentally, these calculations 
reformulate the fact that the Fokker-Planck generators are equal). However, 
the brownian drivers underlying both processes differ, as those for V are 
defined from translations along the N components directions while those for 
V' are associated with the rotation matrices acting on Mr. Interestingly, the 
non-abelian nature of rotation compositions is irrelevant to our discussion. 

Kiessling and Lancellotti [KL06J interpret their representation (A. 5) of 
the Laplace-Beltrami operator on the sphere in terms of two types of con- 
tributions. Here, one associates with three components, 3k — 2<i< 3k, 
the cartesian components of the velocity \'f of particle k in R 3 , so that 
in the Qy's one distinguishes the terms associated with different particles 
([//3] 7^ 17/3] for the ceiling function [•]), and terms coupling two ve- 
locity components of the same particle ( \i/3~\ = \j/3~\ ) so that the particle 
velocity vector just rotates in M? as under a gyroscopic force. As they ob- 
serve, the gyroscopic contribution to each particle velocity drops out in the 
N" -> oo limit. 

In a forthcoming work, we shall return from this model with ,/V" = N' — 1 
effective particles to the more physical model, with particle interactions 
conserving energy and momentum. We shall also eliminate the self-coupling 
gyroscopic terms in Q. and show that, in the limit N" — > °°, one recovers the 
same limit. 
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